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Solar coronal mass ejections (CMEs) show a large variety in their kinematic properties. CMEs originating in active regions 
and accompanied by strong flares are usually faster and accelerated more impulsively than CMEs associated with filament 
eruptions outside active regions and weak flares. It has been proposed more than two decades ago that there are two 
separate types of CMEs, fast (impulsive) CMEs and slow (gradual) CMEs. However, this concept may not be valid, since 
the large data sets acquired in recent years do not show two distinct peaks in the CME velocity distribution and reveal 
that both fast and slow CMEs can be accompanied by both weak and strong flares. We present numerical simulations 
which confirm our earlier analytical result that a flux-rope CME model permits describing fast and slow CMEs in a 
unified manner. We consider a force-free coronal magnetic flux rope embedded in the potential field of model bipolar 
and quadrupolar active regions. The eruption is driven by the torus instability which occurs if the field overlying the flux 
rope decreases sufficiently rapidly with height. The acceleration profile depends on the steepness of this field decrease, 
corresponding to fast CMEs for rapid decrease, as is typical of active regions, and to slow CMEs for gentle decrease, as is 
typical of the quiet Sun. Complex (quadrupolar) active regions lead to the fastest CMEs. 
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1 Introduction 

Coronal mass ejections (CMEs) are spontaneous ejections 
of plasma and magnetic flux from the inner solar corona 
into interplanetary space. They can originate in active re- 
gions, where they are associated with X-ray or EUV flares 
and often also with filament eruptions, but they can also 
originate from eruptions of large quiescent filaments out- 
side active regions not accompanied by strong flares. Early 
observations, based on relatively small data sets, indicated 
that flare-related CMEs are faster, often moving already at 
nearly constant speed (mostly in excess of ^ 750 kms~^) 
when they enter the field of view of white-light corona- 
graphs, i.e., they are fully accelerated at heliocentric dis- 
tances < 2Rq , whereas CMEs with weak or no flare activity 
are accelerated only gradually to moderate speeds (typically 
~ 400-600 km s~^) within the coronagraph field of view of 
6 Rq or more. 

Such data led MacQueen & Fisher (1983) to suggest the 
existence of two separate types of CMEs, flare-associated 
ones and eruption- associated ones (Fig. 1), the latter refer- 
ring to quiescent filament eruptions. Their concept has been 
supported by Sheeley et al. (1999), who categorised the two 
types as impulsive and gradual, respectively. Both author 
groups suggested that the acceleration is dominated by dif- 
ferent physical processes, which may imply the need for dif- 
ferent models. 
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Studies of the large CME data sets meanwhile available 
revealed a very wide range of kinematic properties, with ve- 
locities ranging from ~ 50 - 3000 km s~^ and accelerations 
ranging from ~ 10 to > 10^ ms~^ over ~ 5-1000 min- 
utes. However, the data always showed a continuous distri- 
bution of velocities with a single peak (Zhang & Dere 2006 
and references therein). Furthermore, both fast and slow 
CMEs can be accompanied by both weak and strong flares 
(Vrsnak, Sudar & Ruzdjak 2005). These findings suggest 
that the concept of physically different CME types might 
not be valid, rather there may exist a continuous distribution 
of CME kinematics, with fast impulsive and slow gradual 
CMEs representing the two ends of the distribution. 

Theoretical models of CMEs did not address the wide 
range of CME kinematics for many years. Some modellers 
reproduced the trajectories of individual filament eruptions 
and CMEs in good agreement with the observations (e.g., 
Krall et al. 2001 ; Chen & Krall 2003; Torok & KHem 2004, 
2005; WilHams et al. 2005; Chen et al. 2006). A first attempt 
to model impulsive and gradual CMEs in a unified manner 
was presented by Chen & Krall (2003). These authors en- 
forced the expansion of a line-tied, semi-toroidal flux rope 
by driving it out of equilibrium through the "injection" of 
poloidal magnetic flux at the photospheric level. The accel- 
eration of the rope is dominated by its Lorentz self-force 
(hoop force) and significant only during the period of flux 
injection. A range of acceleration profiles of the rope apex, 
encompassing impulsive and gradual CMEs, was obtained 
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Fig. 1 Left: observation of two ap- 
parently distinct CME types (data from 
MacQueen & Fisher 1983): fast CMEs 
(red) acquire nearly constant speeds low 
in the corona (behind the coronagraph's 
occulting disk indicated in grey), while 
slow CMEs (blue) are accelerated only 
gradually to reach the gravitational es- 
cape speed (dashed line). Right: the TI 
of a freely expanding current ring yields 
similar velocity profiles for different de- 
cay indices n of the external poloidal 
field Bex (equilibrium field strengths at 
the initial radius Ro = 10 Mm are given). 



by prescribing different amplidudes and durations of the in- 
jection and by varying geometrical parameters of the rope. 

Recently, we suggested a flux rope model which ex- 
plains the wide range of CME kinematics without the need 
of such external driving (Kliem & Torok 2006, hereafter 
KT06). As in the model of Chen & Krall (2003), the flux 
rope expands due to the hoop force, but the eruption is 
caused by an ideal MHD instability, the torus instability 
(TI). A range of flux rope trajectories encompassing im- 
pulsive and gradual events results naturally from the model 
if the height dependence of the magnetic field overlying 
the flux rope is varied. However, the analytical treatment 
of CME acceleration in KT06 required simplifications, no- 
tably a freely expanding circular current ring (rather than a 
line-tied flux rope) was considered, and flux pileup, pressure 
gradients, and external toroidal fields had to be neglected. 

In this paper, we present MHD simulations of the TI 
which include the effects of photo spheric line tying and of 
flux pileup in front of the expanding rope. By considering 
a substantial range of magnetic field profiles above the flux 
rope, we confirm our previous result (KT06) that the TI of 
a flux rope provides a CME model which describes both 
impulsive and gradual CMEs and their association with the 
main magnetic structural properties of their source region. 

2 Torus instability 

Observations of erupting filaments, which often evolve into 
the cores of CMEs, suggest the magnetic topology of a 
single arched magnetic flux rope, or partial current ring, 
whose footpoints are anchored in the solar surface. As is 
well known from fusion research, the required pre-eruption 
equilibrium in a low-beta environment must include an ex- 
ternal poloidal magnetic field Bq^, since the hoop force, 
as well as the net pressure gradient force, of a bent non- 
neutralized current channel always point radially outward. 
Such a current ring is unstable against expansion if its hoop 
force decreases more slowly with major ring radius R than 
the opposing Lorentz force due to (Bateman 1978). As- 
suming oc R~^, where n = —RdlnBQ^/dR is the 



decay index of the external field, Bateman (1978) found 
n > Tier = 3/2 as condition for the instability. 

KT06 suggested the instability as a possible initiation 
and driving mechanism for CMEs, referring to it as torus 
instability, and treated its evolution for the first time (using 
the same external field as Bateman 1978). It was found that 
the acceleration profile of the freely expanding current ring 
strongly depends on the steepness of the external field de- 
crease with R. For moderate to steep decrease, n > 2, the 
acceleration rises quickly to a peak within R/Ro ^ 2 and 
decreases quickly afterwards, resulting in a profile as ob- 
served in impulsive CMEs. The peak acceleration increases, 
and its radial position decreases, with increasing n, so that 
the CME develops a more impulsive trajectory if it orig- 
inates in more localized and more complex flux concen- 
trations. For very gradual decrease (n approaching ricr), 
the acceleration has a gradual profile, nearly uniformly dis- 
tributed over a large radial range, i.e., with small peak value 
shifted to larger R; see Fig. 1 in KT06. 

Hence, the instability produces not only a continuum of 
acceleration profiles ranging from fast impulsive to slow 
gradual ones, it also reproduces the observed association 
of impulsive profiles with highly concentrated flux in the 
source (active regions, high n) and of gradual profiles with 
the gradual field decrease over a large height range in the 
quiet Sun (n ~ 3/2, Vrsnak et al. 2002). Fig. 1 shows a 
comparison of velocity profiles obtained from this model 
with observed CME velocities. We note that more gradual 
profiles result also if Rq is increased (which compresses the 
X axis of the right panel in Fig. 1 proportionally); this is the 
effect used by Chen & Krall (2003). Both small n and large 
Ro make large-scale quiescent filament eruptions gradual. 

3 Numerical simulations 

In order to confirm the analytical results, we perform 3D 
zero-beta ideal MHD simulations of the TI, using the ap- 
proximately force-free coronal flux rope model by Titov & 
Demoulin (1999) as initial condition. The line-tied flux rope 
corresponds to a current ring partly submerged in the dense 
subphotospheric plasma. 
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Fig. 2 Initial configurations of the numerical simulations, cor- 
responding to the red, blue, green, and black lines in Figure 3, 
respectively (from top left to bottom right). The rainbow-coloured 
field lines show the centre of the flux rope (current ring), the green 
field lines show the overlying magnetic field (which is a superposi- 
tion of the external poloidal field Bex and the poloidal field created 
by the ring current). Normalized distances between the monopoles 
are 3.6, 2.7, and 0.5 /lo in the bipolar configurations and 2.5 and 
5.0 /io in the quadrupolar one. 



The setup and initial conditions of the simulations are 
similar to those described in Torok & Kliem (2005). Here 
we use a flux rope with a smaller aspect ratio by setting its 
minor radius to 6/10 of its initial apex height ho (to obtain 
a configuration with a flux rope twist below the threshold 
of the helical kink instability) and we choose the initial flux 
rope shape to be almost semi-circular by setting the depth 
of the torus center below the photosphere to be only 0.1 ho 
(since the line-tying suppresses the TI before a semi-circular 
shape is reached; Vrsnak, Ruzdjak & Rompolt 1991; Chen 
et al. 2006). Furthermore, we remove the toroidal compo- 
nent of the external field used in Torok & Kliem (2005), to 
facilitate the comparison with the analytical TI model. The 
resulting purely poloidal external field is created by a pair of 
subphotospheric monopoles (see Titov & Demoulin 1999). 
It is more realistic than the one used in KT06, as it yields 
a decay index which increases monotonically with height 
above the photosphere, n = n{h). 

We fix the geometrical flux rope parameters and con- 
sider different external field profiles by varying the number, 
position and strength of the subphotospheric monopoles in 
the model (such that the variations of the ring current and 
the flux rope twist stay below 1%). The resulting configura- 
tions are shown in Fig. 2. The field above the flux rope of the 
configuration in the bottom right panel is quadrupolar and 
contains an X-line. The top panel in Fig. 3 shows the cor- 
responding decay indices n{h) of the external fields. The 
external field in the quadrupolar configuration drops very 
rapidly with height, with n{h) > 3 already in the region im- 
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Fig. 3 Decay indices n(h) — —hdlnBe^/dh of the external 
poloidal field and corresponding acceleration and velocity profiles 
for the torus-unstable flux ropes shown in Fig. 2 (|n| ^ oo at the 
null in the quadrupolar configuration). The initial position of the 
current-carrying flux rope is indicated in grey. The inner edges of 
the MK3 and LASC0/C2 coronagraph fields of view are indicated 
for a scaled initial apex height of the flux rope of 50 Mm. Veloci- 
ties are normalized to the initial Alfven velocity Va at the flux rope 
apex; Ta = ho/va- The acceleration during the initial relaxation 
phase is omitted for clarity of the figure. 



mediately above the rope. The other configurations all have 
external fields with n{h) < 3 above the rope, but with dif- 
ferent slopes. As ^ oo, n ^ 3 for all configurations. 

The TI is found to occur also in the numerical model. It 
is triggered by moderate upward motions set up in the ini- 
tial relaxation phase of the only approximatively force-free 
configuration, which lift the flux rope to a height where the 
conditions for TI onset (2ho ^ D, D being the footpoint 
distance, and n{h) > ricr) are fulfilled. In the case shown in 
the top left panel in Fig. 2, the field decrease is so gradual 
that the initial relaxation does not move the flux rope to the 
height required for TI onset, the rope finds a stable numeri- 
cal equilibrium after relaxation. In order to trigger the TI in 
this configuration, we artificially set the monopole strength 
to be 9% smaller than its equilibrium value, which causes 
the rope to relax towards an equilibrium at a slightly greater 
height, sufficient for TI onset. The simulations show agree- 
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ment with the analytically obtained instability threshold: in 
the two runs with the lowest n{h) profile, the TI starts at 
flux rope apex heights where n < 1.5. 

Despite the differences between the analytical and the 
numerical model, the rise profiles are qualitatively very sim- 
ilar (compare Fig. 3 with Fig. 1 in KT06). In both cases, the 
acceleration rises quickly to a maximum which strongly in- 
creases with the steepness of the external field decrease. Af- 
terwards, the steeper the field decrease, the steeper the drop 
of the acceleration. In the numerical model, the peak accel- 
eration is shifted to somewhat greater heights for fast ejec- 
tions, h/ho ^ 2-4, which appears to conform better to the 
observations than the analytical result (h/ho ~ 2). 

The quadrupolar configuration yields the most impul- 
sive acceleration and the highest terminal velocity. 

4 Conclusions and discussion 

We present numerical simulations of an erupting force- 
free magnetic flux rope embedded in the potential field of 
model bipolar and quadrupolar active regions. The eruption 
is driven by the TI which occurs if the external poloidal field 
decreases sufficiently rapidly with height. The simulations 
show that the TI occurs not only in freely expanding cur- 
rent rings (KT06), but also in line-tied flux ropes. The TI 
threshold is similar in both cases, ncr ^3/2. 

The acceleration profile depends on the steepness of the 
field decrease, corresponding to fast CMEs for rapid de- 
crease (as is typical of active regions) and to slow CMEs 
for gradual decrease (as is typical of the quiet Sun). These 
simulations confirm our previous finding (KT06) that the TI 
permits to describe fast and slow CMEs in a unfied manner. 
This demonstrates, in accordance with the observational re- 
sults in Vrsnak, Sudar & Ruzdjak (2005) and Zhang & Dere 
(2006), that the concept of two separate CME types is not 
valid and that the wide range of CME kinematics can be 
explained by a single underlying physical process. 

The ratio of the maximum velocities (accelerations) in 
the fastest and the slowest eruption in Fig. 3 is ^ 4 40). 
These ranges are about one order of magnitude smaller than 
the corresponding ranges observed in CMEs. However, it is 
important to note that the peak velocities and accelerations 
listed in the Introduction are matched with very reasonable 
choices for the Alfven speed in the core of large active re- 
gions (Va ~ 3750 kms~^) and for ho (^ Rq/10), although 
the inclusion of line tying and flux pileup in front of the 
expanding rope in the numerical model reduce the accel- 
eration considerably in comparison to the analytical result. 
The great majority of observed CME speeds and acceler- 
ations is accommodated for if the variations of Va and ho 
between 6 spots and quiescent filaments, both by factors 3- 
10, are taken into account. It is always possible to match 
extremely slow CMEs by setting the parameters of the sim- 
ulation very close to the instability threshold. Moreover, for 
gradual CMEs which result from eruptions of quiescent fil- 
aments and are accelerated over several Rq, the presence 



of the heliospheric current sheet reduces the decay index of 
the overlying field in comparison to the asymptotic value 
(n = 3) in the Titov & Demoulin equilibrium. 

Magnetic fields in multipolar active regions decrease 
in the lower corona faster with height than those in bipo- 
lar active regions or in the quiet Sun. Consequently, our 
model predicts the fastest CMEs to originate in multipo- 
lar active regions, which is in line with observations. No 
"breakout-like" reconnection above the flux rope is required 
in our simulation of the quadrupolar model configuration, 
the eruption is solely driven by the TI. 

In contrast to the model of Chen & Krall (2003), our 
model allows to obtain a wide range of CME kinematics 
without imposing an external driving of the configuration 
with a duration similar to the observed duration of signifi- 
cant acceleration. For given flux rope parameters, the kine- 
matics are determined by the slope of the magnetic field 
overlying the flux rope. Although the magnetic field in the 
corona cannot be measured directly at present, a wide range 
of height profiles of the field can be expected to exist, as the 
large variety of field distributions revealed by photospheric 
magnetograms indicates. 

We note that, for a given overlying field, the acceleration 
profile also depends on the chosen flux rope parameters as 
well as on the existence of a flux rope velocity at TI onset 
(due to external perturbations; Schrijver et al. 2007). Fur- 
thermore, we expect that an external toroidal field and pres- 
sure gradients will influence the evolution as well. These 
effects will be studied in forthcoming work. 
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